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1. System for estimating the ground condition under a driving 
vehicle, comprising : 

a wheel speed sensor (4) for sensing a wheel speed signal 
( t(n) , G)(n) ) which is indicative of the wheel speed of a 
vehicle's wheel driving over the ground (2,3); and 
a first analyser unit (8) coupled to said wheel speed sensor 
(4) which comprises 

a sensor imperfection estimation section (9) which is 

designed to estimate a sensor imperfection signal ( S t ) 
from the wheel speed signal ( t(n) ) which is indicative of 
the sensor imperfection of the wheel speed sensor (4) ; 
- a signal correction section (10) which is designed to 
determine an imperfection-corrected sensor signal (s(n)) 
from the wheel speed signal (t n ) and the sensor 
imperfection signal ( 8 X ) ; and 

a ground condition estimation section (11) which is 
designed to estimate a first estimation value ( r(ri) , a(n) ) 
indicative of the ground condition from the imperfection- 
corrected sensor signal ( e(ri) ) . 

2. The system of claim 1, wherein the wheel speed sensor (4) 
comprises a segmented rotary element (5) , and the sensor 
imperfection estimation section (9) is designed to estimate, at 
each revolution of the rotary element (5) , a sensor 
imperfection value ( 8 z ) representative of the sensor 
imperfection signal for each of the segments (6) of the rotary 
element (5) . 

3. The system of claim 2, wherein the sensor imperfection value. 
(<y z ) is a weighted average of sensor imperfection values ( y(n) ) 
of previous and current revolutions (n) of the rotary element. 

4. The system of one of the preceding claims, wherein the 
sensor imperfection estimation section (9) comprises a low pass 
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filter which is implemented according to the following filter 
relation: 

LP: §,=(l-fi)S,+My(n), 

with 

^) = -=^«»)-*(»-l))~ 

wherein is an estimation value of the sensor imperfection, 

ju is a forgetting factor of the filter, r(^) and t(n-\) is the 
wheel speed signal, L is the total number of segments (6) of 
the rotary element (5) and 7^ P (n) is the duration of a complete 
revolution of the rotary element (5) . 

5. The system of one of the preceding claims, wherein the 
ground condition estimation section (9) comprises: 

- a variance determination section (12) which is designed to 
determine the variance (a(n)) of the imperfection-corrected 

sensor signal ( e(ri) ) , and 

- a ground condition estimation subsection (13) which is 
designed to estimate the first estimation value ( r(ri) ) on the 
basis of the variance (a(n)) thus determined. 

6. The system of one of claims 2 to 5, wherein the variance 
determination section (12) comprises a low pass filter (16) for 
determining the variance (a(n)) of the imperfection-corrected 
sensor signal ( s{n) ) according to the following relation: 

a(n) = Var{s) = LP{s 2 ) - LP(s) 2 , 
wherein LP(s) is a low pass filtered value of the imperfection- 
corrected sensor signal ( s(n) ) and LP(s 2 ) is a low pass filtered 
value of the square (s 2 (ri)) of the imperfection-corrected 
sensor signal (s(n)) • 

7. The system of claim 6, wherein the low pass filter (16) is 
implemented according to the following filter relation: 

LP: a{n + \) = (\-X)a(n) + Ae(n), 
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wherein a is an estimation value of the variance Var{s) , X is 
a forgetting factor of the filter, and s{n) is the 
imperfection-corrected sensor signal. 

8. The system of one of the preceding claims, wherein the 
ground condition estimation subsection (13) comprises a signal 
change determination section (14) which is designed to 
determine signal change values ( CUSUMCounter(ri) ) according to 
the following relation: 

CUSUMCounter(n + 1) = min(mzK(CUSUMCounter(n) + a(n) - Drift ,0), CounterLimit) , 
wherein a(n) is the variance obtained from the variance 
determination section, and Drift and CounterLimit are tuning 
parameters . 

9. The system of claim 8, wherein the ground condition 
estimation subsection (13) further comprises a decision section 
(15) which is designed to compare the signal change values 
(CUSUMCounter(n)) from the signal change determination section 
(14) with a first and a second threshold value (set,reset) and to 
output a current first estimation value ( r(ri) ) indicative of a 
rough road condition if the current signal change value 
( CUSUMCounter{ri) ) is greater than the first threshold value 
(set), a current first estimation value indicative of a normal 
road condition if the signal change value ( CUSUMCounter(n) ) is 
lower than the second threshold value ( reset ) , and otherwise a 
current first estimation value equal to the previous first 
estimation value ( r(n - 1) ) . 

10. The system of one of the preceding claims, which comprises: 
- one first analyser unit (8) for each wheel (f = FL,FR,RL,RR ) 

of the vehicle having more than one wheel, wherein each 
first analyser unit (8) is designed to provide a first 
estimation value (a^n)) indicative of the ground condition 
under the respective wheel, and 

a combination section (17) which is designed to combine the 
first estimation values (a z .(n)) provided from each of the 
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first analyser units (8) in order to obtain a combined first 
estimation value ( K w ) ' AiC") ) indicative of the road 
condition under the vehicle. 

11. The system of claim 10, wherein the combined first 
estimation value {y(n) , I hl (n)) is determined by 

- averaging the first estimation values (a.(n)) provided from 
each of the first analyser units (8) , 

- using networks of series expansion type, in particular 
neural networks, radial basis function networks, • fuzzy 
networks, on the basis of the first estimation values 
( a i (n) ) / 

- using a min-function on the basis of the first estimation 
values {aXn)) * and/or 

- using a max-function on the basis of the first estimation 
values (a f (n) ) • 

12. The system of claim 10 or 11, in combination with claim 8 
or 9, wherein the signal change determination section (14) is 
coupled to the combination section (17) in order to determine 
the signal change value ( CUSUMCounter(ri) ) on the basis of the 
combined first estimation value {y(ri)) . 

13. The system of one of the preceding claims, further 
comprising : 

a second analyser unit (19) which is associated with the 
wheel speed sensor (4) and designed to determine a second 
estimation value ( /3(n)) indicative of the ground condition 
from the wheel speed signal {a>{n)) received from the wheel 
speed sensor (4) ; and 

a decision unit (2 0) which is designed to determine a 
combined estimation value (R(n)) indicative of the ground 
condition on the basis of the first and second estimation 
values (a(n) , J3(n)) from the first and second analyser units 
(8,19) , respectively. 
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14. The system of claim 13, wherein the second analyser unit 
(19) comprises: 

a band pass or high pass filter section (21) for filtering 
the wheel speed signal ( a>{n) ) , and a variance estimation 
5 section (12) for determining a variance value (/?(«)) from 

the filtered wheel speed signal ( 3(n) ) , wherein the variance 
value ( fi{n) ) is indicative of the ground condition under the 
respective wheel; 

a side-wise correlation section which is designed to 
10 correlate the wheel speed signals ( co(n) ) of the wheels 

( i = FL,FR,RL,RR ) on a first side of the vehicle (1) with the 
wheel speed signals ( co{n) ) of the wheels (i = FL,FR,RL,RR) on 
a second side of the vehicle (1) , wherein the correlation 
value {r(n)) is indicative of the ground condition; 
is - an axle-wise correlation section which is designed to 
correlate the wheel speed signals ( €o(n) ) of the wheels 
(i = FL,FR,RL,RR) on a first axle of the vehicle (1) with the 
wheel speed signals (co(n)) of the wheels (i = FL,FR,RL,RR) on 
a second axle of the vehicle (1) , wherein the correlation 
20 value (r(n)) is indicative of the ground condition; or 

a frequency determination section which is designed to 
determine the highest Fourier frequency (r(n)) of the wheel 
speed signal ( o)(n) ) which is indicative of the ground 
condition. 

25 

15. The system of claim 13 or 14, comprising: 

- one first analyser unit (8) for each wheel (/ = FL,FR 7 RL,RR ) 
of the vehicle having more than one wheel, wherein each 
first analyser unit (8) is designed to provide a first 
30- estimation value ia^n)) indicative of the ground condition 

under the respective wheel, and 

a first combination section (17) which is designed to 
combine the first estimation values (a^n)) provided from 
each of the first analyser units (8) in order to obtain a 
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combined first estimation value {y(ri)) indicative of the 
road condition under the vehicle; 

a signal change determination section (14) which is designed 
to determine signal change values ( CUSUMCounter(ri) ) on the 
basis of the combined first estimation values (y(ri)) 
according to the following relation: 

CUSUMCounter(n + 1) = mm(max(CUSUMCounter(n) + y(ri)- Drift, 0), CounterLimit) , 
wherein Drift and CounterLimit are tuning parameters; 
one second analyser unit (19) for each wheel (i = FL,FR,RL,RR ) 
of the vehicle, wherein each second analyser unit (19) is 
designed to provide a second estimation value (/?,-(»)) 
indicative of the ground condition under the respective 
wheel , and 

a second combination section (17) which is designed to 
combine the second estimation values (/^(/i)) provided from 
each of the second analyser units (19) in order to obtain a 
combined second estimation, value (r 2 (n)) indicative of the 

road condition under the vehicle 
- an output combination section (22) for combining the signal 
change values (CUSUMCounter(n)) and the second combined 
estimation values (r 2 (n)) in order to obtain a combined 
estimation value ( Q.(n) , R(n) ) indicative of the road 
condition under the vehicle. 

16. The system of claim 13 or 14, comprising: 

one first analyser unit (8) for each wheel (i = FL,FR,RL,RR ) 
of the vehicle having more than one wheel, wherein each 
first analyser unit (8) is designed to provide a first 
estimation value (a^n)) indicative of the ground condition 
under the respective wheel, and 

a first combination section (17) which is designed to 
combine the first estimation values (a f (/?)) provided from 
each of the first analyser units (8) in order to obtain a 
combined first estimation value (r^ri)) indicative of the 
road condition under the vehicle; 
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one second analyser unit (19) for each wheel ( i = FL,FR,RL,RR ) 
of the vehicle, wherein each second analyser unit (19) is 
designed to provide a second estimation value ( /?,(#)) 
indicative of the ground condition under the respective 
wheel, and 

a second combination section (17) which is designed to 
combine the second estimation values ( fi.(n)) provided from 
each of the second analyser units (19) in order to obtain a 
combined second estimation value (r 2 (ri)) indicative of the 
road condition under the vehicle 

an output combination section (22) for combining the first 
and second combined estimation values ( r } (ri) , r 2 (n) ) in order 
to obtain a combined estimation value ( Cl(n) ) indicative of 
the road condition under the vehicle; and 

a signal change determination section (14) which is designed 
to determine signal change values ( CUSUMCounter(ri) ) on the 
basis of the combined estimation values ( Q(ra) ) from the 
output combination section (22) according to the following 
relation: 

CUSUMCounter(n + 1) = min(max(CUSUMCounter(n) + Cl(n) - Drift \0\CounterLimit) , 
wherein Drift and CounterLimit are tuning parameters. 

17. The system of claim 15 or 16, further comprising a decision 
section (15) according to claim 9. 

18. Method for estimating the ground condition under a driving 
vehicle, comprising the steps of: 

sensing a wheel speed signal ( t(n) , co(n) ) by means of a wheel 
speed sensor (4) which is indicative of the wheel speed of a 
vehicle's wheel driving over the ground (2,3); and 

estimating a sensor imperfection signal (S t ) from the wheel 
speed signal ( t(ri) ) which is indicative of the sensor 
imperfection of the wheel speed sensor (4) ; 
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determining an imperfection-corrected sensor signal (a(n)) 
from the wheel speed signal (t(n)) and the sensor 
imperfection signal (<f 7 ); and 

estimating a first estimation value ( r(n) , a(n) ) indicative of 
5 the ground condition from the imperfection-corrected sensor 

signal ( e(ri) ) . 

19. The method of claim 18, wherein the step of estimating the 
sensor imperfection signal from the wheel speed signal 

10 (t(n)) comprises estimating, at each revolution of the rotary- 
element (5), a sensor imperfection value {S t ) representative of 

the sensor imperfection signal for each of the segments (6) of 
a rotary element (5) . 

15 20. The method of claim 19, wherein the sensor imperfection 
value (S t ) is a weighted average of sensor imperfection values 
(y(n)) of previous and current revolutions in) of the rotary 
element . 

20 21. The method of one of the preceding claims, wherein the step 
of estimating the sensor imperfection signal ( 8 } ) from the 
wheel speed signal (t(n)) comprises a step of low pass 
filtering according to the following filter relation: 

LP : S^il-ju^i+Myin) > 

25 wherein 

J<») = -=T^W-^-l))~ 

wherein § l is an estimation value of the sensor imperfection, 
ju is a forgetting factor of the filter, t(n) and t(n-Y) is the 
wheel speed signal, L is the total number of segments (6) of 
30 the rotary element (5) and 7^ F (w) is the duration of a complete 
revolution of the rotary element (5) . 
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22. The method of one of the preceding claims, further 
comprising the steps of: 

- determining a variance (a(n)) of the imperfection-corrected 
sensor signal ( e(ri) ) , and 

- estimating the first estimation value ( r(n) ) on the basis of 
the variance ia(n)) thus determined. 

23. The method of one of claims 19 to 22, wherein the step of 
determining a variance (a(tt)) of the imperfection-corrected 
sensor signal ( e{ri) ) comprises the step of low pass filtering 
the imperfection-corrected sensor signal (s(n)) according to 
the following relation: 

a(n) = Var{s) = LP(s 2 )-LP(s) 2 , 
wherein LP(s) is a low pass filtered value of the imperfection- 
corrected sensor signal ( £■(«)) and LP(s 2 ) is a low pass filtered 

value of the square (s 2 (n)) of the imperfection-corrected 
sensor signal 

24. The method of claim 23, wherein the low pass filtering is 
implemented according to the following filter relation: 

LP : a{n + 1) = (1 - X)a(n) + Xs{ri) , 
wherein a is an estimation value of the variance Var{s) , X is 
a forgetting factor of the filter, and e(n) is the 
imperfect ion- corrected sensor signal. 

25. The method of one of the preceding claims, further 
comprising the step of determining signal change values 
{CUSUMCounter(ri)) according to the following relation: 

CUSUMCounterin + 1) - min(max(CUSUMCounter(n) + cc(n) - Drift, 0), CounterLirnif) , 
wherein a(n) is the variance obtained from the variance 
determination section, and Drift and CounterLimit are tuning 
parameters . 



26. The method of claim 25, further comprising to compare the 
signal change values (CUSUMCounter(n)) with a first and a second 
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threshold value (set, reset) and to output a current first 
estimation value ( r(n) ) indicative of a rough road condition if 
the current signal change value (CUSUMCounter(n)) is greater 
than the first threshold value (set), a current first 
5 estimation value indicative of a normal road condition if the 
signal change value (CUSUMCounter(ri)) is lower than the second 
threshold value (reset), and otherwise a current first 
estimation value equal to the previous first estimation value 
(r{n-l) ) . 

10 

27. The method of one of the preceding claims, further 
comprising : 

- providing a first estimation value («.(«)) indicative of the 
ground condition under the respective wheel for each wheel 
15 (z= FL,FR 7 RL,RR) of the vehicle having more than one wheel, and 

combining the first estimation values in order to 

obtain a combined first estimation value (y(ri) , I hl (n)) indicative 
of the road condition under the vehicle. 

20 28. The method of claim 27, wherein the combined first 
estimation value ( y(n) , I hl (ri) ) is determined by 

averaging the first estimation values (a^n)) provided from 
each of the first analyser units (8) , 

using networks of series expansion type, in particular 
25 neural networks, radial basis function networks, fuzzy 
networks, on the basis of the first estimation values (a^ri)) , 

using a min- function on the basis of the first estimation 
values (a^ri)), and/or 

using a max- function on the basis of the first estimation 
3 0 values (a.(rc) ) . 

29. The method of claim 27 or 28, in combination with claim 8 
or 9, wherein a signal change value ( CUSUMCounter(n) ) is 
determined on the basis of the combined first estimation value 
35 (y(n) ) - 
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30. The method of one of the preceding claims, further 
comprising : 

determine a second estimation value (j3(n)) indicative of the 
ground condition from the wheel speed signal ( co{n) ) received 
from the wheel speed sensor (4) ; and 

determining a combined estimation value (R(n)) indicative of 
the ground condition on the basis of the first and second 
estimation values ( a(n) , /3{n) ) ■ 

31. The method of claim 30, further comprising: 

filtering the wheel speed signal ( a>(n) ) with a band pass or 

high pass filter, and determining a variance value (j3(n)) 
from the filtered wheel speed signal ( 5i(n) ) , wherein the 
variance value ( /3(n) ) is indicative of the ground condition 
under the respective wheel; 

correlating the wheel speed signals ( o)(n) ) of the wheels 
(i = FL,FR,RL,RR) on a first side of the vehicle (1) with the 
wheel speed signals ( a>(n) ) of the wheels (i = FL,FR,RL,RR ) on 
a second side of the vehicle (1) , wherein the correlation 
value (r(n)) is indicative of the ground condition; 
correlating the wheel speed signals ( a>(n) ) of the wheels 
(i = FL,FR,RL,RR) on a first axle of the vehicle (1) with the 
wheel speed signals ( a>{n) ) of the wheels ( i = FL 3 FR,RL,RR ) on 
a second axle of the vehicle (1) , wherein the correlation 
value (r(n)) is indicative of the ground condition; or 
determining the highest Fourier frequency {r(n)) of the wheel 
speed signal ( co(n) ) which is indicative of the ground 
condition. 

32. The method of claim 30 or 31, comprising the steps of: 
providing a first estimation value (a^n)) indicative of the 
ground condition under the respective wheel, for each wheel 
(i = FL,FR,RL,RR) of the vehicle having more than one wheel; 
and 
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combining the first estimation values (ct f (w)) in order to 
obtain a combined first estimation value (/(h)) indicative 
of the road condition under the vehicle; 

determining signal change values ( CUSUA4Counter(n) ) on the 
basis of the combined first estimation values (y(n)) 
according to the following relation: 

CUSUMCounter(n + 1) = min(max(CUSUMCounter(n) + y(ri) - Drift, 0\CounterLimit) , 
wherein Drift and CounterLimit are tuning parameters; 
providing a second estimation value ( j3.(n)) indicative of 
the ground condition under the respective wheel, for each 
wheel ( i — FL 9 FR,RL,RR ) of the vehicle; and 

combining the second estimation values (/?.(/?)) in order to 

obtain a combined second estimation value (r 2 (n)) indicative 
of the road condition under the vehicle; 

combining the signal change values ( CUSUMCounter(ri) ) and the 
second combined estimation values (r 2 (n)) in order to obtain 
a combined estimation value ( Q(«) , R(ri) ) indicative of the 
road condition under the vehicle. 

3. The method of claim 30 or 31, comprising: 

for each wheel (i = FL,FR,RL,RR ) of the vehicle having more 

than one wheel, providing a first estimation value (a^ri)) 
indicative of the ground condition under the respective 
wheel ; and 

combining the first estimation values (a f (n)) in order to 
obtain a combined first estimation value indicative 
of the road condition under the vehicle; 

for each wheel ( i = FL,FR,RL,RR ) of the vehicle, providing a 
second estimation value (jS^n)) indicative of the ground 
condition under the respective wheel; and 

combining the second estimation values in order to 

obtain a combined second estimation value (r 2 (n)) indicative 
of the road condition under the vehicle 
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combining the first and second combined estimation values 
(r x (ri) , r 2 (ri) ) in order to obtain a combined estimation value 
(Cl(n)) indicative of the road condition under the vehicle; 
and 

5 - determining signal change values ( CUSUMCounter(ri) ) on the 
basis, of the combined estimation values ( ) according to 
the following relation: 

CUSUMCounter{n + 1) = min(max(CUSUMCouriter(n) + Q(n) - Drift, 0), CounterLimit) 
, wherein Drift and CounterLimit are tuning parameters . 

10 

34. The method of claims 32 or 33, further comprising the steps 
of claim 26 . 

35. A computer program including program code for carrying out 
15 a method, when executed on a processing system, of estimating 

the ground condition under a driving vehicle, the method 

comprising the steps of: 

sensing a wheel speed signal (t(n) , oo(n)) by means of a wheel 
speed sensor (4) which is indicative of the wheel speed of a 
20 vehicle's wheel driving over the ground (2,3); and 

estimating a sensor imperfection signal {S t ) from the wheel 
speed signal (t{ri)) which is indicative of the sensor 
imperfection of the wheel speed sensor (4) ; 

determining an imperfection-corrected sensor signal (s(n)) 
25 from the wheel speed signal ( t(n) ) and the sensor 

imperfection signal (£;) ; and 

estimating a first estimation value (r(n) , a(n)) indicative of 
the ground condition from the imperfection-corrected sensor 
signal ( s(n) ) - 
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